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ABSTRACT 

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and ^^CO 
(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an 
outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper Hmit for the 
dense core bolometric luminosity of 0.05 Lq is obtained. By modeling the broadband SED and the continuum 
interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the 
observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or 
by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between 
these two models, which produce similar fits. We also detect ^^CO (2-1) emission with red- and blue-shifted 
emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found 
towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the 
best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. 
Whatever the true nature of the central object in LI 451 -mm, this core presents an excellent laboratory to study 
the earliest phases of low-mass star formation. 

Subject headings: ISM: clouds — ISM: individual (L1451, Perseus) — stars: formation — stars: low-mass — 
ISM: molecules 



1. INTRODUCTION 

Star formation takes place in the densest regions of molec- 
ular clouds, usually referred to as dense cores. The parental 
molecular clouds show highly supersonic velocity disper- 
sions, while the dense cores show subsoni c levels of turbu- 
lence ( [Goodman eTaLlp^98l [Caselli et aL][2Q02) ). Recently, 
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Pineda et aT] pO 1 0| ) showed that this transition in velocity dis- 
persion is extremely sharp and it can be observed in NH3 (1,1) 
(see also Pineda et al. 2011, in prep). 

Starless dense cores represe nt the initial conditions of star 
formation. [ Crapsi et"aL] ((2005 ) identify a sample of starless 
cores which show a number of signs indicating that they may 
be "evolved" and thus close to forming a star. 

In the earliest phases of star formation a starless core un- 
dergoes a gravitational collapse. Increasing central densities 
will result in an increase in dust optical depth and thus cooling 
within the core will not be as efficient as in the earliest phases. 
This increases the gas temperature and generates more pres- 
sure. The first numerical simulation to study the formation of 
a protostar from an isothermal core (Larson 1969 ), revealed 
the formation of a central adiabatic core, defined as a "first 
hydrostatic core" (hereafter FHSC). This FHSC would then 
accrete more mass and undergo adiabatic contraction until H2 
is dissociated, at which point it begins a second collapse un- 
til it forms a "second hydrostatic core," which is the starting 
point for protostellar objects. 

A few FHSC cand idates have been suggested in the past. 
Belloche et al.| ( |2006i ) present single dish observations of the 
Cha-MMSl dense core which combined with detections at 24 
and 70 jim with Spitzer suggest the presence of a first hypo- 
static core or an extremely young protostar (see also Bell ochej 
|et aL|2011| ). |Chenet al.| ( |2010 ) present SMA observationToF 
the continuum at 1.3-mm and ^^CO (2-1) line in the L1448 
region located in the Perseus cloud where no Spitzer (IRAC 
or MIPS) source is detected. They detect a weak contin- 
uum source and a wel l collimated high- v elocity outflow is ob 
served in i^CO (2-1). 



Chen et al. ( 2010 ) analyze different sce- 
narios to explain the observations and conclude that a FHSC 
provide the best case, however, no actual modeli ng of the in 



terfero metric observations is presented. Recently, [Enoch et aL 
|2010|) present CARMA 3-mm continuum and deep Spitzer 



70 jim observations of another FHSC candidate (Per-Bolo 58) 
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in the NGC1333 region also located in the Perseus cloud. In 
these observations t hey detect a weak so urce in the 3-mm con- 
tinuum and 70 /im. [Enoch etH] pOl Q) simultaneously mod- 
eled the broadband SED and the visibilities, allowing them to 
conclude that the best explanation for the central source is a 
FHSC. Dunham et al. (2011, submitted) present SMA 1.3- 
mm observations which reveal a collimated slow molecular 
outflow using ^^CO (2-1) emission. 

Another class of low luminosity objects has been identified 
thank s to Spitzer: Very Low Luminosity Objects (VeLLOs 
e.g. Young et al ||20Q4| [Bourke et aT][2n05t IDunham et al. 



some of which are found within evolved cores (as clas 



sified by |Crapsi et al.|20Q5] ). These objects have low intrinsic 
luminosities (L < 0.1 Lq) and are embedded in a dense core 
( |di Francesco et al.||2007|). As VeLLOs have only recently 
been revealed by Spitzer dPunham et al.| 2008 ), it is not yet 
clear whether these are sub- stellar objects that are still form- 
ing, or low-mass protostars in a low-accretion state. 

Broadband SED modeling of VeLLOs suggest that these 
sources can be explained as embedded YSOs with a surround- 
ing disk. In the case of IRAM 04191+1522 (hereafter IRAM 
04191), continuum observations using the IRAM Plateau de 



Bure in terferometer (PdBI) were interpreted by Belloche et al. 
( |2002| ) as produced from the dense core's inner part without 
the need for a disk. 



Recently, |Maury et al.| ( [2010| ) presented high-resolution 
PdBI observations towards a sample of 5 Class sources to 
study the binary fraction in the early stages of star forma- 
tion. Their sample includes two previously known VeLLOs: 
L1521-F and IRAM 04191. Dust continuum emission is de- 
tected toward both objects, which may arise from either a cir- 
cumstellar disk or from the inner parts of the envelope. Lack 
of detailed modeling of the SED or visibilities in these sources 
makes it hard to distinguish between these two scenarios. 

This paper presents observations of LI 451 -mm, a low-mass 
core without any associated mid-infrared source in which we 
have detected compact thermal dust emission and a molecular 
outflow, along with models constructed to derive the proper- 
ties of this object. In ^we discuss previous observations of 
LI 451 -mm. In ^presents data used in this paper. In ^we 
present the analysis of the observations and radiative transfer 
models to reproduce the observed spectral energy distribution 
(SED) and the continuum visibilities to constrain the physical 
conditions of the source. Finally, we present our conclusions 
infl 

2. L1451-mm 



L1451-mm (also known as Per-Bolo 2; "Eno ch et al.|2006) ) 
is a cold dense core in the LI 451 Dark Cloud located in the 
Perseus Molecular Cloud Complex. Here we assume that 
Perseus is at a distance of 250 pc ( Cerni sjl990t |Hirota et al.| 
|2008), which is consistent with those used by previous works. 
LI 451 -mm is detected in 1.1 mm dust continuum with Bolo- 
cam at 31'^ resoluti on, and its estimated mass is 0.36 Mq from 
the Gaussian fit by [Enoch et al.| ( [2Q06| ) with major and minor 
FWHMs of 33^' and 54'', respectively. However, the core is 
too faint to be identifi ed by the SCUBA surveys at 850 /xm 
of the Perseus Clou d ( jHatcheU et al.|[2Q05l |Kirk et al.||2006l 
[Sadavoy et al.|2010|). 

Figure [TJpresents a summary o f the observations pre - dating 



this worOowards L1451-mm. [Foster & Goodman] ( |2006 
presented deep Near-IR observations (J H Ks) of L1451-mm 
which show only heavy obscuration, and no evidence for a 
point source. Establishing upper limits for this non-detection 



Table 1 

Photometry of L 1 45 1 -mm 



Filter 


Wavelength 
(12m) 


Flux 
(mJy) 


Aperture 
(arcsec) 


/ 


1.25 


< 0.006 


2.25 


H 


1.65 


< 0.0048 


2.25 


Ks 


2.17 


< 0.0099 


2.25 


IRACl 


3.6 


< 0.048 


2.5 


IRAC2 


4.5 


< 0.012 


2.5 


IRAC3 


5.8 


< 0.060 


2.5 


IRAC4 


8.0 


< 0.030 


2.5 


MIPSl 


24.0 


< 1.5 


16 


MIPS2 


70.0 


<72 


16 


MIPS3 


160. 


880 ±330 


16 


IRAM 


1200 


70±7 


16.8 


Note. — 


Upper limits used are 3-cr limits. 



was complicated by the presence of extended bright structure 
(i.e., cloudshine) around the edge of L1451-mm. We esti- 
mate an upper limit by inserting synthetic stars with a range of 
magnitudes (in 0.1 magnitude steps) and appropriate FWHM 
at the central position. We ran Source Extractor ( Berlin & 
[Arnouts 1996 ) on these synthetic images using a 2.25'' radius 
aperture and established the input magnitude at which sl 3 a 
source was successfully extracted. 



This core is classified as "starless" by [Enoch et al.| 

because no poin t source is detected in Spitzer IRAC and MIPS 
images (J0rgens en et al.| 2006: Reb ull et al.|2007|). Since the 
IRAC images do not contain significant extended emission 
we measured the flux in a 2.5^' radius aperture centered on 
the central position of LI 451 -mm with a background annu- 
lus of 2.5 to 7.5^' using the IRAF phot routine and applied 
the aperture correction factor for this configuration from the 
IRAC instrument handbook. All fluxes measured this way 
were within 2a of zero (fluxes were both positive and nega- 
tive). For MIPS we used the smallest aperture with a well- 
defined aperture correction factor, which is 16'^ Both MIPSl 
and MIPS2 were consistent with zero flux while MIPS 3 was a 
weak (2.1a) detection. A summary of the photometric results 
is presented in TablelT] 

Given the lack of detectable emission at Spitzer wave- 
lengths, and using the correlatio n between 70/im and i ntrinsic 
YSO luminosity determined by [Dunham et al.| ( [2008| ), an up- 
per limit of L < 1.6 X 1O~^L0 on the luminosity of a source 
embedded within LI 451 -mm is determined. 

For a given SED, two quantities can be calculated to de- 
scribe it: bolometric luminosity, L^^/, and bolometric temper- 
ature, 7^0/ • The bolometric luminosity is calculated through 
integration of the SED {Sj,) over the observed frequency 
range. 



i-bol =^7Td I Sjydu 



(1) 



while the bolometr ic temperature is calculated following |My- 
|ers & Ladd| ( [T993] ), 



where 



r^^/ = 1.25x10"^^ 



Hz 



(2) 



(3) 



For LI 451 -mm, if the upper limits are used as measurem ents, 
then we obtain Ljjoi < 0.05 and Tiyoi < 30K (see ,Dunham| 
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Figure 1. Summary of observations available towards L1451-mm between 1 and 160/im. (a) Ks, H and / three color image, red, green, and blue respecti vely. 
The white cross shows the central position observed with the SMA. (b) lRAC-4, -2 and -1 three color image, red, green, and blue respectively, (c) M lPI^T 
monochrome, (c) MlPS-2 monochrome; and (d) MlPS-3 monochrome hgures. No point source is detected in any ot the images. 



let al.t 2Q0F; Enoch et al.|2009b| for discussions on the uncer- 
tainties in calculating T^oi andL/,^/). This bolometric lu minos- 
ity is lowe r than any of the Class objects studied by |Enoch 
et al.| ( [2QQ9b| ) in Serpens, Ophiuchus and Perseus Molecular 



Table 2 

VLA Spectral Setup: 23 GHz Setting 



Clouds; and also it is fainter than any of the VeL LQs with 
(sub-)niillimeter wavelength observations studied by Dunham 
|t al. (2008). 

Pineda et al. (2011 in prep) present NH3 (1,1) and (2,2) 
line maps observed with the 100-meter Green Bank Tele- 
scope. From these observations they derive an almost constant 
(within a ^ 1' radius) kinetic temperature, T^in ^9.1 K, and 
velocity dispersion, ^ 0.15kms~\ showing no evidence 
for heating from a central source. 

3. OBSERVATIONS 

3.1. Single dish continuum observations 

Dust continuum observations at 1.2-mm were taken us- 
ing MAMBO at IRAM 30m telescope, under good weather 
(^i.2mm = 0.1-0.2). The data reduction was carried out us- 
ing MOPSI, with parameters optimized for extended sources. 
The observations are convolved with a 15^' Gaussian kernel, 
while the flux unit is in Jy per 1 1^' beam. The rms noise level 
is 1 mJy per 11^' beam, and the map for the core studied is 
shown in Figure [2] 

3.2. VLA Observations 

Observations were carried out with the Very Large Array 
(VLA) of the National Radio Astronomy Observatory on Jan- 
uary 10, 2006 (project AA300). The NH3 (J,K) = (1,1) and 









Resolution 


Frequency 


Molecule 


Transition 


Chan. 


(kms-i) 


(GHz) 


NH3 


(1,1) 


63 


0.1544 


23.694495 


NH3 


(2,2) 


63 


0.1542 


23.722733 



(2, 2) inversion transitions were observed simultaneously (see 
Table |2] for a summary of the correlator configuration used). 
At this frequency the primary beam of the antennas is about 
1.9^ The array was in the compact (D) configuration, the 
bandwidth was 1.56 MHz, and the channel separation was 
12.2 kHz (corresponding to 0.154 kms~^). This configuration 
is centered at the main hyperfine component and it also covers 
the inner pair of satellite lines for NH3 (1,1). 

The bandpass and absolute flux calibrator was the quasar 
0319+415 (3C84) with a calculated flux density of 10.6 Jy at 
1 cm, and the phase and amplitude calibrator was 0336+323. 
The raw-data were reduced using CASA image processing 
software. The signal from each baseline was inspected, and 
baselines showing spurious data were removed prior to imag- 
ing. The images were created using multi-scale clean (scales 
[8,24,72] arcsec and smallscalebias=0.8) with a robust param- 
eter of 0.5 and tapering the image with a 8^' Gaussian to in- 
crease the signal-to-noise. Each channel was cleaned sepa- 
rately according to the spatial distribution of the emission, us- 
ing a circular beam of 8'^ Table |3] lists relevant information 



4 



J. E. Pineda et al. 



on the maps used. 

3.3. CARMA observations 

Continuum observations in the 3 -mm window were ob- 
tained with CARMA, a 15 -element interferometer consist- 
ing of nine 6.1 -meter antennas and six 10.4-meter antennas, 
between April and September 2008. The CARMA correla- 
tor records signals in three separate bands, each with an up- 
per and lower sideband. We configured one band for maxi- 
mum bandwidth (468 MHz with 15 channels) to observe con- 
tinuum emission, providing a total continuum bandwidth of 
936 MHz. The remaining two bands were configured for 
maximum spectral resolution (1.92 MHz per band) to observe 
NH2D (lii-loi) and N2H+ (1-0) (see Table|4]for the corre- 
lator configuration summary). The six main hyperfine com- 
ponents of NH2D fit in the two narrow spectral bands and 
six of the seven-hyperfine components of N2H"^ (1-0) were 
observed, with the highest frequency (isolated) component 
falling outside the observed frequency range. 

The field of view (half-power beam width) of the 10.4- 
m antennas is 66^' at the observed frequencies. Seven 
point mosaics were made around the center of L 1451 -mm 
in CARMA's D and E-array configurations, giving baselines 
that range from 8-m to 150-m. Observations of N2H"^ (but 
not NH2D) were also made in CARMA's C-array configura- 
tion, with projected baselines of 30-m to 350-m. The syn- 
thesized beam sizes and position angles (measured counter 
clockwise from North) are: 5.4^' x 4. 8^' and -77° (continuum), 
5.2''x4.3'' and -73° (N2H+), 9.2''x7.6'' and 72.3° (NH2D). 
The largest angular size to which these observations were sen- 
sitive is ^40". 

The observing sequence for the CARMA observations was 
to integrate on a primary and secondary phase calibrator 
(3C 111 and 0336+323) for 3 minutes each and the science 
target for 14 minutes. In each set of observations 3C 111 
was used for passband calibration and observations of Uranus 
were used for absolute flux calibration. Based on the repeata- 
bility of the quasar fluxes, the estimated random uncertainty 
in the measured fluxes is a 5%. Radio pointing was done at 
the beginning of each track and pointing constants were up- 
dated at least every two hour s thereafter, us ing either radio or 
optical pointing routines (Corder et al.'2010). Calibration and 
imaging w ere done us ing the MIRIAD data reduction pack- 
age ( Sault et al.|1995) . Table |3] lists relevant information on 
the maps used. 

3.4. SMA observations 

The SMA observations were carried out at 1.3 -mm 
(230 GHz) in both compact and extended configuration. The 
compact array observations were carried out on Novem- 
ber 1, 2009, with zenith opacity at 225 GHz of -0.085. 
Quasars 3C 84 and 3C 1 1 1 were observed for gain calibra- 
tion. Flux calibration was done with observations of Uranus 
and Ganymede. Bandpass calibration was done using obser- 
vations of the quasar 3C 273. The SMA correlator covers 
2 GHz bandwidth in each of the two sidebands. Each band 
is divided into 24 "chunks" of 104 MHz width, which can be 
covered by varying spectral resolution. The correlator config- 
uration is summarized in Table [5] 

The extended array observations were carried out on 
September 13, 2010, with zenith opacity at 225 GHz of 
—0.05. Quasars 3C 84 and 3C 111 were observed for gain 
calibration. Flux calibration was done with observations of 



Uranus and Callisto. Bandpass calibration was done using ob- 
servations of the quasar 3C 454.3. The SMA correlator used 
the new 4 GHz bandwidth in each of the two sidebands. Each 
band is divided into 48 "chunks" of 104 MHz width, which 
can be covered by varying spectral resolution. The correlator 
configuration is summarized in Table |6] 

Both data sets were edited and calibrated using the MIR 
software packagq^ adapted for the SMA. I maging was per- 
formed with the MIRIAD package dSault e t al. 1995 ), result- 
ing in an angular resolution of 1.23 xO. 88'^ R\=85.6° (using 
robust weighting parameter of -2) and 1 .35'' x 0.96'' PA=80.9° 
(using robust weighting parameter of 0) for the continuum and 
^^CO (2-1), respectively. Table[3]lists relevant information on 
the maps used. The rms sensitivity is ^0.5 mJy beam~^ for the 
continuum, using both sidebands (avoiding the chunk contain- 
ing the ^^CO line), and — 36mJybeam~^ per channel for the 
line ^^CO (2-1) data. The primary beam FWHM of the SMA 
at these frequencies is about 55". 

4. RESULTS 
4.1. MAMBO and CARMA continuum 

The MAMBO dust continuum emission map (left panel of 
Figure [2]) can be decomposed into a bright compact core, and 
fainter filamentary emission. The compact core peak position 
is located at {a, 6) =(03:25:10.4, +30:23:56.0). The compact 
core within 4,200 AU mass is estimated to be 0.3 M^:), where 
a dust opacity per dust mass of 1.14cm^g~^ (jOssenkopf &| 
|HenningH1994| ), gas-to-dust ratio of 100, and a dust temper- 
ature of 10 K are used. The MAMBO derived mass is con- 
sistent with the mass previously estimated using Bolocam. In 
Figure [3] the compac t core is compared to the sample of star- 
less cores from Kauf fmann et al.| ( |2008| ), where the fiducial 
radius of 4, 200 AU is used to compare with previous works 
(e.g., Motte & Andre 20 01). suggesting that it is more com- 
pact than most starless cores From this comparison we can es- 
tablish that LI 451 -mm is in fact quite compact, and therefore 
dense, suggesting that it is more comp act than most starl ess 
cores an evolved evolutionary state (see Crapsi et al.|20"Q5] ). 

A faint central source is detected in the CARMA 3-nmi con- 
tinuum map shown in right panel of Figure |2] The contin- 
uum emission map is fitted by a Gaussian witn a total flux of 
lOmJy, while if a point source is fitted a flux of (4±2) mJy is 
obtained. A summary of the fits to the CARMA 3 -mm contin- 
uum are listed in Tabled The CARMA continuum emission 
agrees with the MAMBO peak position. 

4.2. Molecular lines with VLA and CARMA 

Figure |4] shows the summary of the molecular line transi- 
tions observed with CARMA and VLA. In this study we will 
briefly discuss the kinematics of the region and leave a more 
in depth study of the core in forthcoming papers (Schnee et 
al., 2011 in prep., Arce et al., 2011 in prep.). From these ob- 
servations a centroid velocity and ve locity dispersion are ob- 
tained by fitting the line profiles (see |Rosolowsky et al.|2008| 
for details). The integrated intensity maps show the extended 
emission from the core where the peaks match the position 
of the CARMA continuum emission to within the respective 
beam size. The centroid velocity maps, for all three lines, 
show a consistent result with a clear velocity gradient. A gra- 
dient is fitted to the centroid velocity map for all three lines, 

|http : //cf a- www . harvard. edu/$\sim$c qi7| 
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Table 3 

Parameters of Interferometric Maps 



Map Array Beam^ rms 



NH3 (1,1) 


VLA 


8''x8''(0°) 


3mJybeam~^ channel"^ 


NH3 (2,2) 


VLA 


8''x8''(0°) 


3mJybeam~^ channel"^ 


NH2D(ln-loi) 


CARMA 


9.2''x7.6''(-h72.3°) 


90mJybeam~^ channel"^ 


N2H+ (1-0) 


CARMA 


5.2''x4.3''(-73°) 


80mJybeam~^ channel"^ 


3 -mm continuum 


CARMA 


5.4^^x4.8'' (-77°) 


0.5 mJy beam" ^ 


i^CO (2-1) 


SMA 


1.35''x0.96''(-h80.9°) 


40mJybeam~^ channel"^ 


1.3 -mm continuum 


SMA 


1.23''x0.88''(-h85.6°) 


0.5mJybeam~^ 



^ Size and position angle. Position angle is measured counter clockwise from north. 




Figure 2. Left panel shows the MAMBO 1.2 mm dust continuum emission, where there is both a compact central bright object (Speak = 33 mJybeam ^) and also 
less bright and diffuse emission. Black solid contours represent [1, 2, 3, ... , 15] x 2.5 mJybeam"^ levels, and dashed contours are -2.5 and -5 mJybeam"^ levels. 
The dashed rectangle shows the area imaged by CARMA data. Right panel shows the gray scale 3-mm continuum emission map observed with CARMA. Solid 
contours mark the [2, 3, 5] xO.5 mJybeam"^ while negative contours are shown by dotted lines. A faint central source is detected at the center of the image that 
matches the MAMBO peak position. 



Table 4 

CARMA Spectral Setup: 3-mm Setting 



Table 6 

SMA (extended configuration) Spectral Setup: 1.3-mm Setting 









Resolution 


Frequency 


Molecule 


Transition 


Sideband Chan. 


(kms-i) 


(GHz) 


NH2D 


lii-loi 


Lower 2 x 63 


0.106 


85.9262 


N2H+ 


1-0 


Upper 2x63 


0.098 


93.1737 






Table 5 






SMA (compact configuration) Spectral Setup: 1.3-mm Setting 








Resolution 


Frequency 


Molecule 


Transition 


Chunk Chan. 


(kms-i) 


(GHz) 




LSB 


ci^o 


2-1 


s23 512 


0.28 


219.560357 


13CO 


2-1 


sl3 256 


0.55 


220.398684 


USB 


12CO 


2-1 


sl4 256 


0.53 


230.537964 


N2D+ 


3-2 


s23 512 


0.26 


231.321966 



Note. — For all other chunks the channels have a resolution of 
0.8125 MHz. 



with an average value of Q =6.1 kms~^ pc~^ and a -66° posi- 
tion angle (measured counter clockwise from north), see Ta- 
ble [8] for the individual fit obtained for all three maps. This 



Molecule 


Transition 


Chunk Chan. 


Resolution 
(kms-i) 


Frequency 
(GHz) 




LSB 


Ci8o 

13CO 


2-1 
2-1 


s23 128 
sl3 128 


1.11 
1.11 


219.560357 
220.398684 


USB 


12CO 

N2D+ 


2- 1 

3- 2 


sl4 256 
s23 128 


0.53 
1.05 


230.537964 
231.321966 



Note. — The channels of all other chunks have resolution 0.8125 MHz, 
except those in chunks sl5 and sl6 where the resolution is 1.625 MHz. 



velocity gradient is larger than those observed in low er angu- 
lar resolution NH3 (1,1 ) maps ([Goodman et al.|1993| ) or using 
lower density tracers ( [Kirk et al.||2010| ), while velocity gra- 
dients of a similar magnitude are obt ained with high angular 
resolution observation s of dense gas jCurtis &" Richer 20^ 
Tanner & Arce"2011). The velocity dispersion maps show a 
clear increase towards the center of the map, starting with very 
narrow lines (close to the thermal values) in the outer regions. 
The increase in velocity dispersion is more pronounced in the 
N2H"^ velocity dispersion map, and the difference can be ex- 
plained by the higher angular resolution obtained in the N2H"^ 
observations. 
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Table 7 

Results of Fits to CARMA 3-mm continuum image for L1451-mm 





Center^ 


Peak Flux 


Size (FWHM) 


PA 


Source 


a(J2000) (5(J2000) 


(mJy) 


(arcsec) 


(deg) 


Point Source 


3:25:10.25 +30:23:55.09 


4±2 






Gaussian 


3:25:10.21 +30:23:55.20 


2±1 


(16±13, 8±6) 


-37±44 



^ Units of R.A. are hours, minutes, and seconds. Units of declination are degrees, ar- 
cminutes, and arcseconds. 
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••■L1451-mm 
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\0 



;^ 0.1 1.0 

Mass within 4200 AU [Mq] 

Figure 3. Relation between the intrinsic radius at 70% peak intensity and the 
mass within 4200 AU radius from the peak for L1451-mm (show by the solid 
star), compared to the sample of starless cores and VeLLOs candidates from 
[Kauffmann et al. ( 2008). Both properties (mass and radius) are derived from 
MAMBO data. Starless cores with well established and uncertain properties 
are shown by filled and open circles, respectively. Cores hosting candidate 
VeLLOs are shown by stars, the radius bias due to internal heating by the cen- 
tral object is indicated by the arrows (see Kauffmann et al. 2008 for details). 
Curves of constant H2 central density are shown by dotted lines. The dashed 
line indic ates the upper radius limit for evolved dense cores, <4800 AU, sug- 
gested by |Crapsi etaL] |2005). 



Table 8 

Velocity Gradients Fits 



Transition 


Q 


PA 






(km s~^ pc~^ ) 


(deg) 


(kms-i) 


N2H+ 


5.6±0.2 


-84 ±3 


3.970 ±0.002 


NH2D 


8±1 


-83 ±15 


3.949 ±0.009 


NH3 


6.24 ±0.06 


-65.8 ±0.7 


4.0042 ±0.0006 



The NH3 (1,1) and (2,2) integrated intensity maps obtained 
using the VLA are shown in Figure [5] left and middle pan- 
els respectively, where all components observed are taken 
into account. Both lines present a peak coincident with the 
CARMA continuum peak, and where the NH3 (1,1) emission 
covers a more extended region than the NH3 (2,2). However, 
since the emission NH3 (1,1) is fairly extended, the addition of 
GBT data to provide the zero- spacing is needed to allow a ro- 
bust temperature determination. The morphology of the NH3 
and N2H"^ integrated intensity maps show non-flattened struc- 
tures, which are dras tically d ifferent from those seen in young 
Class sources (e.g., [Wi seman et al. 2001 ; Chiang e t al.|2010t 
Tanner &Arce 2011| ). The right panel of Figure [SlpresentstEe 
derived kinetic temperature obtained from the simultaneous 
NH3 (1,1) and (2,2) line fit, with uncertainties in the tempera- 
ture determination between 0.2 K, in the central region, up to 
0.5 K in the outer regions. Surprisingly, the kinetic tempera- 
ture map is quite constant, in particular, there is no evidence 



for an increase in temperature towards the peak continuum 
position. 

4.3. SMA 

The visibility amplitude as a function of wv-distance for the 
1.3-mm continuum is shown in Fig. [6] From this figure we 
identify two components. An extenaed component, that is 
resolved at long baselines, and a compact component that re- 
mains unresolved even at the longest baselines (indicated by 
the horizontal line in Figure [6]). This unresolved component 
is commonly seen towards dense cores containing a central 
protostar, and it is interpreted as arising from an unresolved 
central disk (e.g., J0rgen sen et al.; 2007, 2009). However, in 
the case of LI 451 -mm there is no infrared detection of a cen- 
tral protostar. 

The SMA 1.3-mm continuum map is shown in Figure [Tj 
with red- and blue-shifted ^^CO (2-1) emission overlaid. The 
dust continuum emission clearly shows the central source, 
also presented in Figure [6] The position of this continuum 
source coincides with the pixel where the CARMA 3-mm 
continuum peaks. The orientation of the SMA continuum 
emission, obtained through a fit of the visibilities and listed 
in Table [9j is close to the right ascension axis and clearly 
different from the red- and blue-shifted ^^CO emission. The 
detected ^^CO (2-1) emission shows a large velocity disper- 
sion, with spatially separated blue- and red- shifted lobes (see 
Figure [7]). Figure shows the Position Velocity (PV) dia- 
gram along the grayline drawn in Figure |7] where the dashed 
vertical line shows the centroid velocity of the dense core, 
3.94kms~\ which is consistent with the interferometric ob- 
servations (see Table [8]) and the NH3 (1,1) data obtained with 
the GBT at 30'' (Pineda et al. 2011, in prep). 

The dust mass of the compact (unresolved) emission is es- 
timated as: 



M 



1.3 mm ■ 



13 Mr, 







1.3 mm 



200 pc 



)(' 



,0.36(30 K/T)_ 



')■ 

(4) 

where we assumed optically thin emission and a dust opacity 
per dust mass (a^i j mm) of 0.86 cm^ (thick ice man tles co- 
agulated at 10^ cm~^ from'Ossenk opf & Henning|l994i ) and a 
gas-to-dust ratio of 100, see J0rgen sen et al.l ( |2007| ). The flux 
from the unresolved emission is estimated by fitting a p oint 
source to baselines longer than 40 ^A, as in J0rgensen et al.| 
( |2007| ), and therefore it avoids contamination from the dense 
core itself. The result of fitting a point source gives a flux of 
27.0 mJy, see Table [9j which implies a mass of 

(5) 



M 



1.3 mm 



= 0.024 Mc7 







for a temperature of 30 K (as used in J0rgensen et al.||2007| ). 
This mass will be us ed as a fir st estimate for the circumstellar 
disk mass (see J0rgensen et al.,2007^ for a discussion). 
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03'25'^12' 11^ 10' 9' 

RA (J2000) 




RA (J2000) 

Figure 4. Top, middle, and bottom rows present results from N2H+ (CARMA), NH2D (CARMA), and NH3 (VLA) emission line maps, respectively, where all ob- 
served hyperfine components are used. Left, middle, and right columns show the integrated intensity, centroid velocity, and velocity dispersion maps, respectively. 
Left panels also show contours for the integrated intensity at the following levels: [5,15,25] Jybeam~^kms~^ for N2H+ (top), [5,15,25,35,45,55] Jybeam~^kms~^ 
for NH2D (middle), [12.5,50,87.5,125,162.5,200] mJybeam~^kms~^ for NH3 (bottom). The color scales for the velocity (centroid and dispersion) maps are in 
kms~^ . Contours show the CARMA 3-mm continuum emission presented in Figure[2] The synthesized beam for each transition line is shown at the bottom left 
corner and the CARMA 3-mm continuum beam shown above. The orange box in the upper left panel shows the region imaged using the SMA. 



Table 9 

Results of uv Fits of SMA 1.3 -mm continuum data for LI 451 -mm 





Phase Center^ 


Offset 


Peak Flux 


Size (FWHM) 


PA 


Source 


a(J2000) (5(J2000) 


(arcsec) 


(mJy) 


(arcsec) 


(deg) 


Point Source (vis. longer than 40kX) 


3:25:10.21 -^30:23:55.3 


(0.40, -0.23) 


27.0±0.4 






Gaussian (all vis.) 




(0.40, -0.23) 


32.8±0.6 


(0.66±0.05, 0.45±0.05) 


-88±8 



Units of R.A. are hours, minutes, and seconds. Units of declination are degrees, arcminutes, and arcseconds. 
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(mJy beam"' km s"*) (mJy beam"* km s"*) (K) 



-75 25 125 225 -4.5 0.0 4.5 9.0 7.5 8.5 9.5 




03^25^^12^ 11^ 10' 9^ 
RA (J2000) 



Figure 5. Left and middle panels show the NH3 (1,1) and (2,2) integrated intensity VLA maps of the same region shown in Figure [4] using all the observed 
hyperfine components. The overlaid contours show the integrated intensity at [2,4,6,. . . ,12] x 18.12 and 1.09 mJy beam" ^ km s~^ for NH3 (1,1) and (2,2) in the 
left and middle panels, respectively. The right panel shows the kinetic temperature derived by fitting simultaneously both NH3 lines, overlaid with the CARMA 
3-mm continuum. The kinetic temperature uncertainties ranges between 0.2 K in the central region up to 0.5 K in the outer region. The kinetic temperature map 
presents small variations, where there is no increase in temperature at the peak position of the continuum emission. 
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SMA 1 .3-mm continuum 



: 0.04 P I L1451-r 
0.03 



No signal 
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Figure 6. Visibilities amplitude as a function of wv-distance. The dotted 
histogram indicates the expected amplitude in the absence of signal. The 
solid line shows the flux from the point source fit reported in Table|9] 
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Figure 7. The greyscale map shows the source detected in the 1.3-mm con- 
tinuum continuum observed with SMA, with the overlaid black contours at 
[3,7,11,15,19,23] mJybeam-i. Overlaid are the contours for the ^^CO (2-1) 
integrated intensity using the red and blue channels (red and blue channels 
are taken between 5.3-6.9 km s~^ and 1.9-3.7 km s~^ respectively). Dotted 
contours denote negative contour levels. Contour levels are drawn at inte- 
ger multiples of 114mJybeam~^ kms~^. Dashed light grey contours show 
the CARMA 3-mm continuum emission presented in Figure |2] The white 
line shows the direction of the Gaussian fit on the uv-plane as reported in 
Tablej9| The gray line is cut for the position-velocity diagram shown in Fig- 
urelSPrhe 1.3-mm continuum ^^CO (2-1) emission synthesized beams are 
shown at bottom left and right corners, respectively. 



The disk-to-dense core mass ratio, M disk /M dense core^ is es- 
timated using the disk mass from eq. [5] and the dense core 
mass. This ratio is low (^0.1) but comparable to Class ob- 
jects dEnoch et al.||2Qllt |j0rgensen et al.||2Q09[ [Enoch et al. 
|2Q09a| ). 



4.4. Simultaneous Fit of Visibilities and Broadband SED 

A powerful way to constrain the physical parameters of 
dense cores and YSQs is by fitting the broadband SED (e.g., 
Robitaille et al.|[2QQ7] ). In the case of LI 451 -mm, only de- 
tections at 160 and 1200 fim are available, which makes the 
broadband SED fit a not well-constrained problem. Here, the 
information from the 1.3-mm continuum observations (SMA) 
is extremely important to help discri minate between different 
physical models (e.g., Enoch et al.|2009aj ). 

In order to compare the continuum emission model with 



o 
o 
w 
o 



CD 
W 

o 



< 




Figure 8. Position velocity diagram for the SMA ^^CO (2-1) data along 
the gray line shown in Figure [t] Contours drawn start at 0.3 Jybeam"^ with 
an increment of 0.24 Jy beam . The dashed vertical line shows the centroid 
velocity of the dense core, 3.94km s~^. 

the interferometric observations, the model is sampled in uv- 
space to match the observations using the uvmodel task 
in MIRIAD. The synthetic and observed visibilities are both 
binned in uv-distance (using the uvamp task in MIRIAD), 
and then they are added as an extra term to the to mini- 
mize, 

2 _ [ ^i.obs ~ ^i.model \ 



(6) 



where Vi^obs and Vi^modei are the average observed and synthetic 
visibilities in the / bin, respectively, while a^^^^ is the uncer- 
tainty of the observed average visibility. The subject to 
minimization is 



x' = (i-e)xLD+exL. 



(7) 



where Q is an ad-hoc weight used to control how important 
it is to fit the visibilities compared to the broadband SED. In 
this case a 2 = 0.5 is used, which gives the same weight to the 
SED and visibilities fit. 

Because of the problem's high dimensionality, th e min- 
imization is carried out with a genetic algorithm (Oohnston 
et al. 2011), while the model SEDs and visibilities are cal 
culated with a new Monte-Carlo radiation transfer code (Ro- 
bitaille et al., 2011 in prep ), which is based on th e radiative 
transfer code presented by [Whitney et al.| ( [2003b| ). The new 
code uses raytracing for the thermal emission at sub-mm and 
mm wavelengths, providing excellent signal-to-noise to fit the 
long- wavelength SED and visibilities. 

Using this fitting program we explore three models with 
increasing levels of complexity to explain our observations: 
a) starless isothermal dense core; b) dense core with a YSO 
and disk at the center; and c) dense core with a central FHSC. 

The parameter ranges searched using the genetic algorithm 
is given in Table [T0| a summary of the fit results is shown in 
Figure [9j and the model parameters are listed in Table [TT 



Starless isothermal dense core — The simplest model consist 
of a pure isothermal dense core with a density profile nocr'P, 
where the exponent, p, of the density profile is a free param- 
eter, but constrained to values smaller than 2. The dense core 
temperature, T^ust, and the outer radius. Rem, are also parame- 
ters in the fit, which are listed in Table [TT] along with the dense 
core mass, Menv The best model fit forme starless isothermal 
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Starless Isothermal Dense Core 
Dense Core with Disk and YSO 
Dense Core with a central FHSC 
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Figure 9. Summary of best-fit of the broadband SED and dust continuum 
visibilities for three different models: starless isothermal dense core in red; 
dense core with a YSO and disk at the center in blue; and dense core with a 
central FHSC in green. The data are shown in black. The top panel shows the 
visibilities in filled circles. The bottom panel shows the broadband SED for 
L1451-mm, where the upper limits are shown by triangles, measurements are 
shown by filled black circles, and the best model fits are shown by the solid 
curves. 



Table 10 

Parameter Ranges Searched 



Parameter 



Description 



Value/Range Sampling 



Starless Dense Core 



p Exponent of density profile -2-0 linear 

Mdense core Envelope Mass {Mq) 10~^ - 10 logarithmic 

T^ust Envelope temperature (K) 5-30 linear 

Rmax Envelope outer radius ( AU) 1 , 000 - 5 , 000 logarithmic 



Dense Core with Disk and YSO 



Lint 


Intrinsic luminosity {Lq) 


10-4 


-0.1 


lo^ 


^arithmic 


Minfall 


Infall rate(M0yr-i) 


10-9 


- 10-6 


lo^ 


^arithmic 


Rem 


Outer envelope radius (AU) 


1,000 


-5,000 


lo^ 


^arithmic 


^cent 


centrifugal radius (AU) 


50- 


1,000 


lo^ 


^arithmic 


Rdisk 


outer disk radius (AU) 


50- 


1,000 


lo^ 


^arithmic 


^disk 


Disk Mass {Mq) 


10-8 


-10-3 


lo^ 


^arithmic 


i 


Viewing Angle (deg) 


0- 


-90 




linear 



dense core model is shown in red in Figure |9] where it shows 
that this model does not provide a good match to the visibil- 
ities. If the power-law density exponent is not constrained, 
then a very steep density profile, n ex r~^-^, can actually match 
both the SED and visibilities. However, a power-law exponent 
of -2.8 is beyond the range deemed physically reasonable. 
Because, even though the outer region of starless cores (and 



cylin drical filaments) can have similar steep density profiles 
(e.g., |Tafalla et al.||2502| ), in the inner region (r <3,000 AU) 
their density profiles are flat, which is exactly the region we 
are interested in to produce the compact emission. 

Dense Core with a central YSO and disk — The next model fitted 
is one composed of a dense core with a central YSO and disk. 
The dense core is modelled as a rotating and infalling enve- 
lope ( [Ulrich||1976| ), with outer radius 7?env, total mass Mgnv, 
and infall rate Minfaii. The density of the envelope is given in 
spherical polar coordinates by 



p{r,e,4>) = pT 



-3/2 



1 + 



Mo 



-1/2 



ii. 
Mo 



(8) 

where is the centrifugal radius, ji = cos 0, and /io is the 
cosine of the polar angle of a streamline of infalling particles 
as r — > oo, which is given by: 



Mo + Mo 



-M 



:0 



(9) 



The normalization constant p^^^ is related to the infall rate by: 

Mnfall 



^env _ 

Po - 



47r {GM^Rl 



1/2 



(10) 



where is the mass of the central object. The disk is mod- 
eled as a passive flared disk described in cylindrical polar co- 
ordinates by 



(3-q 



exp 



h(R) 



(11) 



where pj^^ is defined by the disk mass Mdisk, Q is the surface 
density radial exponent (which we set to -1), f3 is the disk 
flaring power (set to 1.25), and the disk scale-height h(R) is 
given by: 



h(R) = ho 



R 



lOOAU 



(12) 



where ho is the scale-height at 100 AU and it is set to 10 AU. 

The viewing angle / is a parameter in the fitting. The cen- 
tral protostar is modeled as an object with an effective surface 
temperature of 3,000 K and intrinsic luminosity L^t, where 
the parameter L/„^ includes the luminosity due to accretion. 

The temperature is computed self-consistently with t he den- 
sity using the radiation transfer code, see Whitn ey et al. 



2603b|a |2004J ). It assumes a geometry (e.g., Ulrich envelope 
model with a flared disk), the dust properties, and local t her- 
modynamic equilibri um. Here we use dust opacities of Os-| 
|senkopf & Henning (1 994| ) for dust grains with thick ice man- 
tles after 10^ years of coagulation at a density of 10^ cm~^. 

The best model parameters are listed in Table [TT] This 
model provides an excellent fit to the visibilities, while the 
SED fit underestimates the flux at 160/im. 

Given the best fit result we estimate the accretion luminos- 
ity. Lace, as 



Lacr — 



R^ 



(13) 



where and R^ are the protostellar mass and radius, and 
Mace is the accretion rate onto the protostar. Using equa- 
tion ([13]) we obtain an accretion luminosity of 6Lq, where 
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Table 11 

Best Model Fit Parameters 



Model 


Lint 


Minfall^ 


Mdisk 


Rdisk 


Viewing Angle 


^dense core 


Renv 


Tdust 




(Lq) 


(Moyr-i) 


(Mo) 


(AU) 


(deg) 


(Mq) 


(AU) 


(K) 


Starless Dense Core 












0.353 


1001 


10.2 


Dense Core with Disk and YSO 


0.0450 


7.1 X 10-6 


0.086 


107 


64.9 


0.146 


1007 




Dense Core with central FHSC 












0.56 


2543 





^ Infall rate derived using eq. {14} for M* = 0.086 . 



we have assumed that the accretion rate is the same as the 
infall rate, 



Macr. = M, 



infall ■ 



: 2.4 X 10" 



0.5 



Meyr-\ (14) 



the minimum mass of the central object is the mass of the 
disk, = Mdisk = O.O86M0, and that the central object might 
be a young protostar, = 3Rq. This expected accretion lu- 
minosity is 100 times higher than what can be kept hidden 
at the center of the core in the best-fit model, and therefore 
some of the assumptions must be clearly misrepresenting re- 
ality. Another simple estimate that is derived from equa- 
tion ([13]) is the accretion rate onto the central object needed 
to proHuce the same luminosity of the best model, obtaining 
Mace = 2.45 X lO~^M0yr~^ (much lower than the estimated 
infall rate, 7.1 x lO~^M0yr~^). Clearly, in order to make a 
self-consistent model an extremely low accretion rate must be 
assumed. 

Dense Core with a central FHSC — Another possibility to rec- 
oncile the low luminosity observed and the accretion rates ex- 
pected from the models presented above is to increase the stel- 
lar radius, . There is one object which has been predicted 
from numerical simulations JLarson] 11969] [Masunaga et al. 

^ ' "20001 IMachidc " ' 



T998lpasimaga & Inutsuka|2000[ |Machida et al.|2008| ) that 
is large enough to tit the description : the first hydrostatic core 
(FHSC). 

Masuna ga et al.| fl998) performed radiation hydrodynamic 
simulations of a collapsing core until the formation of a FHSC 
and found that for an initial core of 0.3 Mq (model M3a), a 
FHSC of R^ ~ 5 AU and Lint ^ 0.03 (including accretion 
luminosity) is formed. We use the density and temperature 
profile resulting from this simulation as the input parameter 
for the radiative transfer code. The mass and radius of the 
dense core are varied (see Table [TT]) to find the best match 
of the observed SED and visibilities, see green points in Fig- 
ure |9] Notice that this model does not use the same density 
and temperature profile as used in the dense core with disk 
and YSO model, and this explains the different SEDs. It is 
clear that the FHSC model visibilities do not provide as good 
a match to the observations as the YSO model, however, the 
FHSC model only have two free parameters compared to the 
seven free parameters of the YSO model. 

4.5. The nature of the ^^CO (2-1) emission 

If there is no source of heating within the core, then ^^CO 
should freeze-out onto dust grains (e.g., Tafalla et al. 2002). 
Therefore, the presence of the ^^CO (2-1) emission in itself 
is strongly suggestive of a central heating source. From the 
NH3, N2H"^ and NH2D centroid velocity maps it is clear that 
the velocity gradient is in the RA direction while the ^^CO 
emission is more or less perpendicular. Since the orientation 



of the dust continuum emission detected with SMA is almost 
perpendicular to the ^^CO (2-1) emission, we argue in favor 
of a central object and outflow system. 

The amount of mass needed to keep material at 560 AU with 
a velocity of 1.3 kms~^ (similar parameters to the ^^CO emis- 
sion) is ^ 0.53 M0, which is almost twice the total mass in 
the dense core. Therefore, despite the low velocity seen in the 
^^CO (2-1) emission this gas is unbound and it is consistent 
with ^^CO (2-1) tracing a slow molecular outflow. 

4.6. Outflow Properties 

The physical parameters of molecular outflows a re typically 
calculated using both ^^CO and ^^CO lines (e.g., Arce et al. 
2010). Unfortunately, in L1451-mm there are no detections of 
i^CO (2-1) and we are left to use only the ^^CO (2-1) emis- 
sion to study the high- velocity gas. 

A lower limit for the mass entrained by the outflow, Mfiow, 
is estimated assuming that the ^^CO (2-1) emission is opti- 
cally thin (see details in Appendix |A]) and with an excitation 
temperature of Tex = 20 K. The momentum (P/iow) and energy 
(Efio^) of the outflow along the line of sight are estimated fol- 
lowing [Cabrit & Bertbut, (1990) , 

Pflow = Mput^i I Vi - V center \ , (15) 
i 

Eflow=^ ^Mout4 \Vi-Vcenter\^ , (16) 



where Mout^i is the mass in voxel /, Vcenter is the velocity of the 
core, and v/ is the velocity of voxel /. Also, the outflow charac- 
teristic velocity, Vfiow, is calculated as Pfiow/Mfiow. An upper 
limit for the lobe size, Riobe^ is estimated from the (not decon- 
volved) extension of the red- and blue-shifted ^^CO emission 
(2Riobe) ^ 4.5'', or ^ 1 , 120 AU at the distance of Perseus. We 
also calculate the dynamical time, Tdyn = Riobe/^fiow, mechan- 
ical luminosity, Lfioy, = Efioy,/rdyn, force, Ffioy, = P/iow/rdyn, 

and rate, Mfiow = M flow/ Tdyn- 

The outflow properties are calculated using only the vox- 
els with signal-to-noise ratio higher than 2. To avoid con- 
tamination in the outflow parameters from the cloud emission 
the central channel is not used in the calculations, and a sec- 
ond estimate is calculated where the three central channels are 
removed. The outflow parameters are reported in Table [12 
where we notice that the differences between the quantities 
calculated using both methods are small. 

The outflow properties presented in Table [12] show an ex- 
tremely weak outflow in LI 451 -mm. Howeve r, when com- 
pared to the outflow found towards LI 014 by |Bourke et al. 
( 2005| ) both Pflow and Efiow are close to the lower limits esti- 
mated using similar assumptions. 

The dynamical ti me is cons istent with the FHSC estimated 
lifetime (e.g., .Machida et al.|2008j . 
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Table 12 

Summary of Outflow Properties^ 



Property 



Mass (M0) 

Momentum (Mq km s~^ ) 
Energy (ergs) 
Luminosity (Lq) 
Force (M0kms~^ yr~^) 
Characteristic Velocity (kms~ 
Dynamical Time (yr) 
Outflow rate (Mq yr~^ ) 



1.2 X 10-^ (8.4 X 10-6) 
1.7 X 10-5 (1.4 X 10-5) 
3.1 X 10^^ (3.0 X 10^8) 
1.3x10-6(1.6x10-6) 

8.3 X 10-9 (9.3 X 10-9) 

1.3(1.7) 
2.0 X 10^ (1.6 X 10^) 
6 X 10-9 (5 X 10-9) 



Note. — Properties are calculated without the central 
channel, while the values in parentheses are calculated with- 
out the 3 central channels. 

^ Not corrected for outflow inclination with respect to the 
plane of the sky. 

Recent three-dimensional radiation magneto hydrodynam- 
ics (R MHD) simulatio ns of the de nse co re collapse ( Machidal 

" -"J 



eraLl[2008[ [Tomida et~aLl|201Qbt [Commergon et al.||201()l 
show that when the FHSC is formed a slow outflow can 
be driven even before the existence of a p rotostar (see also 
Tomisaka]|2002[ [Banerjee & Pudritz||20Q6| ). In these simula- 
tions, the outflow that is generated is poorly collimated and 
typically has maximum velocity of 3kms~^- very similar to 
the observed outflow in LI 4 51 -mm. In contrast, theoretical 
models (Shang et al. 2007; Pudritz et al.||2 007 ) and observa- 
tions (e.g., Arce et al. 20 07 ) indicate that outflows from young 
protostars are highly collimated and exhibit velocities of a few 
tens of kms~^ although outflows fro m VeLLOs d isplay lower 
velocities (And re et al. 1999; Bourke et al.|2005| ). 

The properties of the LI 451 -mm molecular outflow are 
consistent with a picture where a first hydrostatic core is the 
driving source of the poorly collimated and slow outflow ob- 
served. 

5. DISCUSSION 

The detection of an unresolved source of continuum in the 
CARMA and SMA observations strongly suggests the pres- 
ence of a central source of radiation and/or a disk. Moreover, 
the simultaneous fit of the broadband SED and the continuum 
visibilities rules out the possibility of explaining the obser- 
vations without a central source (either a YSO or a FHSC). 
From the SED modeling it appears feasible to hide the central 
YSO even at 70/im, but it also requires an extremely ineffi- 
cient or episodic accretion process. This might be consistent 
with the results obtained by Enoch et al. (2009b); Dunham 
et al.| ( |2010) and Offner & McKee (20 11 ), where episodic ac- 
cretion IS argued to explain the low luminosity of YSOs ob- 
served by Spitzer. 

If the unresolved emission observed with the SM A is inter- 
preted as a disk (as done by J0rgens en et al.||2007| ), then the 
disk mass is already ~10% of the dense core mass, w hich is 
similar to the disk mass found in class objects (e.g., [Enoch 



et al.|201 l||2QQ9at|Maury et al.|201Q| ), although jBelloche etal 
(|2002|) shows evidence for a sniall disk in the young protostai 



young protostar 
IRAM 04191 (Mdisk < IO-'Mq and Mcore ^ I.SMq). These 
studies of Class objects suggests that the assembly of mass 
to form a disk starts very early on. 

A slow molecular outflow is detected in the ^^CO (2-1) 
line, see Sec. 4.5 and |4.6| Its orientation is almost perpen- 
dicular to the velocity gradient seen in dense gas tracers ob- 
served with VLA and CARMA (NH3, N2H+ and NH2D), and 



despite the low velocity the gas is unbound. The properties 
presented in Table [T2]place it as the weakest outflow found so 
far, with the lowest energy and momentum measured. Unfor- 
tunately, we have no estimate of the outflow inclination angle 
and therefore some of the outflow parameters might be un- 
derestimated. If the outflow is close to the plane of the sky, 
then the outflow velocity would be faster but it still could be 
consistent with a slow outflow driven by a FHSC depending 
on how large is the correction. This, however, would imply 
that the outflow extension is the one measured in the data, and 
therefore the outflow would have a shorter dynamical time 
and low degree of collimation. On the other hand, if the out- 
flow is nearly in the line-of-sight, then the outflow velocity 
is similar to the 1.5 kms~^ measured from the data. However, 
the outflow extension would be much larger implying a longer 
dynamical tim e, which might be similar to those predicted by 
[Tomida eFaI.| ( |2010a ) for FHSCs in recent numerical simu- 
lations. Therefore, constraining the inclination angle (e.g., 
through observations of the outflow cavity as in [Huard et aL 
2006| ) would provide important insight regarding the outflow 
and by extension to the central object. 

For all the reasons listed above, we claim that a central 
source of radiation (either a YSO or a FHSC) must be present 
within L 1451 -mm. The lack of sensitive observations at mid- 
infrared wavelengths restricts our ability to carry out a more 
detailed modeling of this object. From our best-fit models we 
predict that LI 451 -mm should be detected by the Herschel 
Gould Belt Survey ( Andre & Saraceno 2005), similar to the 
observations by Linz et al . (2010). Therefore, those observa- 
tions will provide a definitive answer regarding the luminosity 
of the central source and give more constraints to the model- 
ing. 

One way to explain our observations is by having a first 
hydrostatic core at the center of the dense L 1451 -mm core, 
instead of a YSO. The simultaneous fit of both visibilities and 
broadband SED shows that a FHSC can also provide a good fit 
to the observations, with the advantage of having an accretion 
luminosity consistent with the observations. The presence of 
a slow and poorly collimated outflow further supports this sce- 
nario. It is for these reasons that we propose L1451-mm to be 
a FHSC candidate. 

Future observations of LI 451 -mm with interferometers us- 
ing a more extended configuration and/or different frequen- 
cies will probe the currently unresolved continuum emission. 
We expect that such observations will provide a constraint on 
the origin of the emission (i.e., disk or first hydrostatic core). 
And, if the disk is confirmed, then a comparison with more 
evolved disks can be carried out. Moreover, observations of 
^^CO (3-2) would provide an estimate of the gas temperature, 
and therefore a good test to confirm that the ^^CO emission is 
generated by an outflow (where the gas is usually warm). 

It is very important to note that three out of the four known 
FHSC candidates are found in the same molecular cloud 
( [Enoch et al.|2010[|Chen et al.J2010| ). We compare this num- 
ber to the expected number of FHSC in Perseus assuming a 
constant star formation rate, which can be estimated as 



FHSC 



Class - 



tfhsc 

' "TClass 



(17) 



We estimate the FHSC lifetime to be ~ 10^ yr (e.g., |Machida 
[et al.[[20Q8]), the number of Class sources in P erseus is 
20-35 ( [Hatchefl et al.[|2007l [Enoch et al. ''2009b). and the 



Class lifetime is 2-5 xWyr ^sser et al.,200L; Hatctiell| 
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et al.|2007 Enoch et al.|2Q09b| ). Finally, the expected number 
of FHSC in Perseus is < 0.2 objects (similar re sults are ob- 
tained using statistics for Class I objects, e.g., |Evans et"aL 
|2b09), and therefore, if all three candidates are confirmed, 
either Perseus is in an extremely peculiar epoch (e.g., a re- 
cent burst on the star formation rate) or this stage is longer 
than previously predicted by numerical simulations. A longer 
lifetime for t he FHSC s tage, up to 10"^ yrs, has recently been 
suggested by |Tomidaet al. (2010a) for FHSCs formed in low- 
mass dense cores (~ 0.1 M©). 

6. SUMMARY 

We present IRAM 30-m, CARMA, VLA, and SMA obser- 
vations of the isolated low-mass dense core L1451-mm in the 
Perseus Molecular Cloud. No point source is detected towards 
the center of the core in NIR and Spitzer observations; how- 
ever, a dust continuum source is identified in both CARMA 
and SMA continuum maps. Upper limits on the bolomet- 
ric luminosity and temperature, Li^oi and Ti^ou of 0.05 and 
30 K are estimated. Also, ^^CO (2-1) emission is observed to- 
wards LI 451 -mm suggestive of a slow and poorly collimated 
outflow. Modeling the broadband SED and observed visibil- 
ities at 1.3-mm confirms the need for a YSO or a First Hy- 
drostatic Core (FHSC) to explain the observations. However, 



more high-resolution observations are needed to distinguish 
between these two scenarios. 

Although YSO and FHSC models are almost indistinguish- 
able, the FHSC scenario seems more likely from the data at 
hand (and thus we may call LI 451 -mm a FHSC candidate). 

Finally, if all current FHSC candidates are confirmed, then 
an important revision of the FHSC lifetime must be carried 
out, which may include modificat ions of the numerical simu- 
lations (e.g., |Tomida et al.|2010a ). 
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APPENDIX 

CALCULATION OF CO COLUMN DENSITY 

If the levels of the molecule are populated following a Boltzmann distribution of temperature Tex, then the column density can 
be expressed as. 



^^y^ gj 1 



Jrdv 



gj+i A 



(1- 



-hv/kTex 



Nj = 93.2H 



2/+1 



Jrdv 



2J + 3J Aj^i^j (l-e-To/Te. 



(Al) 



where gj = (2/+ 1) is the statistical weight of level / for a linear rotor molecule, v is the transition frequency in units of GHz, 
Aj+i^j is the spontaneous emission coefficient in s~^ r is the transition optical depth, the velocity is in kms~\ and 



To = ^^ (A2) 
k 

In the case of ^^CO (2-1), we use u = 230.538GHz and A2^i = 6.91 x 10~^s~^ (obtained from Leiden Atomic and Molecular 
Databas^, and therefore Tq = 11.0641 K. 

Using the equation of radiative transfer to relate the observed emission, Tr, with excitation temperature, T^x, and background 
temperature, Tcmb^ we obtain: 



Tr = To 
idv = To 



1 



1 



(l-O 



j TRdv = To (e^'^^- ~ ^ ~ (^e^o/Tc^' ~ ^ / 



rdv . 



(A3) 



where optically thin emission is assumed. 
Combining equations (|A1|) and (|A3|), the column density of the level / can be calculated as. 



A^/ = 93.28 



2/+1 



/ TMBdv 



M(^^CO) = 9.924x 10^"^ 



2/+ 3/ Aj+l^J (^i_^-To/T,,^ 7q (^e'^o/Te, - 1) ^ - {^e'^o/Tcmb - 1^ 

/ TMBdv 



cm 



(A4) 



where equation ( A4 ) gives the column density of the level / = 1 of ^^CO using the ^^CO (2-1) transition emission. 



http : / / www . strw . leidenuniv . nl/ $\sim$moldata/ 
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The total column density of ^^CO is calculated as 

NC^CO) = ^^NjC^CO) = ^^^e'^'^'^^'/'^-NjC^CO) , (A5) 

gje~^J/'^^^^ (2/+1) 

where B is the rotation constant for a linear rotor (B = 57.635968 GHz for ^^CO), and Z is the partition function, which can be 
approximated as Z kTex/QiB). Therefore, in the case of ^^CO / = 1 we obtain, 

N^CO) = 1.2 X 10- (A6) 

(^ii.06/r,._ 1^-1 _(^^ii.06/r,^,_ 1^-1 (l-^-ii-06/r,.) 

which when combined with a ^^CO abundance with respect to H2, [^^C0/H2] = 10""^, provides an estimate of the total column 
density of H2. 

The final conversion between column density and mass is done using 

^ im) (^) (^)' {^) (^^)"- 

where Asky is the area on the sky used to calculate A'^(^^CO), and /i is the mean molecular weight. 
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